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ABSTRACT 

This article explores the agreement between the predictions of Modified Newtonian Dynam¬ 
ics (MOND) and the rotation curves and stellar velocity dispersion profiles measured by the 
DiskMass Survey. A bulge-disk decomposition was made for each of the thirty published 
galaxies, and a MOND Poisson solver was used to simultaneously compute, from the bary- 
onic mass distributions, model rotation curves and vertical velocity dispersion profiles, which 
were compared to the measured values. The two main free parameters, the stellar disk’s 
mass-to-light ratio (M/L) and its exponential scale-height (hz), were estimated by Markov 
Chain Monte Carlo modelling. The average best-fit K-band stellar mass-to-light ratio was 
M/L ~ 0.55 ± 0.15. However, to match the DiskMass Survey data, the vertical scale-heights 
would have to be in the range hz = 200 to 400 pc which is a factor of two lower than those 
derived from observations of edge-on galaxies with a similar scale-length. The reason is that 
modified gravity versions of MOND characteristically require a larger M/L to fit the rotation 
curve in the absence of dark matter and therefore predict a stronger vertical gravitational field 
than Newtonian models. It was found that changing the MOND acceleration parameter, the 
shape of the velocity dispersion ellipsoid, the adopted vertical distribution of stars, as well as 
the galaxy inclination, within any realistic range, all had little impact on these results. 
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1 INTRODUCTION 


Understanding the dynamics of disk galaxies is essential 
to the vett ing process of theories of galaxy formation and 
ology jFlores & Primack 1^9j^ dTRfok&McGauglJ^ 1^9^ 


de Blok et al.ll200ll:lvan den Bosch & Swatersll200ll:ISwaters et al.l 


20031: iGentile et al.l |2004 lOilmore et all l2007l: Ide Blold l201(]h . 


Disk galaxies moderately inclined to the line of sight (50-80°) 
can provide an HI rotation curve from which the dark matter 
(DM) content and distri bution can be deduced in the context 
of New t onian dynamics ( Bosmal floTSl : iRubin et'^ll978l : iBosmal 
Il981al f9: ISofue & Rubinll200lh . However, there exist degeneracies 


between the DM halo density profile, the stellar mass-to-light ra¬ 
tios of t he stellar components and the scale-height of the stel¬ 
lar disk jvan Alb ada et all fl 9851 : iKuiiken & Gilmore|[l98^ Il99ll : 
lAngus et al.ll2012h. 

The vertical stellar distribution is assumed to be a declining 
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exponential on both sides of the mid-plane. Thus, the scale-height 
of the stellar disk is simply the exponential rate at which the stel¬ 
lar luminosity density drops, with increasing height above or below 
the mid-plane of the disk. It is assumed to be constant with radius, 
which is supported by observations (e .g. Ivan der Kruit & Searld 
ll98lLll982l:fBizvaev & Mitronovall2002h . 

The stellar mass-to-light ratio (M/L) is the ratio between the 
mass of a stellar population and the luminosity, as observed through 
a particular bandpass. It cannot be observed directly, so it is ordinar¬ 
ily inferred by numerical stellar population synthesis models. This 
modelling crucially depends on the star-formation and chemical en- 
richment history considered and on the initial mass function (IMF; 
iBell & de Jong|200ll:IChabriej2003l:l&ouDj200lh . The IMF is the 
spectrum of stellar types formed given a molecular cloud of certain 
initial mass, metallicity and other relevant properties. It therefore 
has the scope to vary from galaxy to galaxy. In addition to the IMF, 
stellar population synthesis models have other source s of uncer¬ 
tainty l lConrov et aI]|2009L I 2 OIOI: IConrov & GuiriilbOlOh . 

In standard fitting of rotation curves of highly inclined disk 
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galax ies, it is common to invoke the “maximum disk hypothe- 
sis”, ( Ivan Albada & Sancisilll98^ : ISacket^ll997l : [Courteau & Rixl 
Il999li i.e. that the stellar disk contributes maximally to the 
rotation curve. This hypothesis is supported by observations 
that deduce the microlensing optical depth in the Milky Way 
(e.g..lBissantz & Gerhard| l2002h. the baryonic Tully-Fisher relation 
dMcGaugh & Schomberlll2015h and measurements that place the 
co-ro tation radius of barred galaxie s just beyond the end of the bar 
(e.g.. ISellwood & Debattistall2014ll . It yields values for the M/L 
that are typically in accordance with the predictions of stellar pop¬ 
ulation synthesis models. 

This, ho wever, does not demonstrate that the hypothesis is cor¬ 
rect (see e.g. iHerrmann & Ciardulldl200^ : iDutton et alJl20TIh and 
it would be ideal to have a robust, independent measurement that 
breaks the disk mass degeneracy. This is theoretically possible be¬ 
cause the M/L can be determined dynamically through the verti¬ 
cal velocity dispersions of stars over the full projected area of the 
galactic disk, if the scale-height is known ( lBahcalJl984ll . For close 
to edge-on disk galaxies, however, one cannot measure the vertical 
velocity dispersion and thus the technique is limited to only those 
galaxies with moderate inclinations to the line of sight (5 — 45°) , 
where 25 < i < 35° is seen as optimal dBershadv et al .11201 Obi : 
hereafter DMSii). 

In order to break the degenerac y between the DM ha lo and 
M/L, the DiskMass Survey (DMS: Is^ershadv et alj|20103 : here¬ 
after DMSi) made observations of the line of sight velocity dis¬ 
persion profiles of 46 nearly face-on disk galaxies (DMSi): 30 of 
which have been published and a further 100 are part of the larger 
survey. They also measured their surface brightness profiles and 
rotation curves. There were two further, essential ingredients in 
the analysis that come from scaling relations. The first is the lu¬ 
minous Tully-Fisher (TF) relation between the absolute magnitude 
of a galaxy and a measure of its outer rotation speed. This helps 
to isolate the inclination of the galaxy, which for nearly face-on 
galaxies can otherwise be ob tained using the tilted disk method of 
lAndersen & Bershadvl l l2013l) . 

The second is the correlation between the disk stellar scale- 
height, a measure of the thickness of the disk, and scale-length, 
a measure of the radial extent. This relationship is explored in 
de tail by DMS i i (the ir §2.2), and is deri ved from the studie s 


. This 


of Kregel et al] 1 2002|) ( hereaft er K02) and iPohlen et al 
ISchwarzkog f&_Dettn^ l l200(]ti : lxilouris et al.l l ll997Lll99' 
allowed Bershadv et all l l201lh and Martinsson et al.l ( I2013al l5l 
(hereafter DMSvi and DMSvii) to infer the M/L of the disk. 
The data imply that the stellar di sks are “sub-maxim al” (K- 
band M/L ~ 0.3 or lower, see ISwaters et all [20141) which 
means they do not contribute maximally to the rotation curve 
in the central regions, contrary to the value found from popula- 
tion synthesis models that as sume a Rroupa IMF (M/L ~ 0.6; 
iMcGaugh & Schombertll2(B^ . This leaves more room for DM in 
the central regions. 

Modified Newtonian D ynamics (MOND rMilgrom|[l983L see 
iFamaev & McGau^l2012l for a recent review) is a theory which 
proposes a modification of dynamics whose impact is most appar¬ 
ent in regions of low acceleration. Most current working versions 
of MOND consist of an actual modification of gravity, i.e., at the 
classical lev el, a modificatio n of the Newtonian Poisson equation 
(but see also lMilgromll20lTI) . This modification occurs due to the 
hypothesised existence of a new constant of physics with dimen¬ 
sions of acceleration, ao ~ For accelerations much 

stronger than this threshold, ao, there is no discerned deviation 
from Newtonian gravity. However, far below the threshold the true 


acceleration perceived by a test mass is found from — gwao - 
whe re qn is the expecte d Ne wtonian gravitational fi eld. 

iNipoti et ^ 1 I 2 OO 7 I) and iBienavme et al.l l l2009li made studies 
of the vertical dynamics of the Galaxy in MOND, showing that 
it could be possible to distinguish between MOND and the DM 
paradigm with data from the Milky Way. The extra constraint on 
the dynamics from vertical velocity dispersions in nearly face-on 
disks of external galaxies provides a new test of this hypothesis. 
This article addresses whether MOND can simultaneously account 
for the measured vertical velocity dispersions and rotation curves, 
while keeping in line with galaxy scaling relations. 

In section 2 the framework is presented for the joint modelling 
of galaxy rotation curves and stellar vertical velocity dispersions in 
the MOND context. In section 3 the methods are discussed and this 
includes a discussion of the bulge-disk decomposition, the accu¬ 
racy of the Poisson solver, and the observational error budget for 
the main data. In section 4 the primary results are presented, this 
includes a discussion of the fits to the vertical dynamics and rota¬ 
tion curves, the confidence ranges of the fitted parameters, and how 
well the fitted parameters mesh with other observations. In section 
5 possible scenarios that could alter the results are discussed. In 
section 6 conclusions are drawn and their implications are explored. 


2 DYNAMICAL ANALYSIS OF THE DISKMASS SURVEY 
2.1 Rotation curve fitting 

The following reviews how to fit the measured rotation curve 
of a disk galaxy in an idealised case. This is done in order to expose 
the free parameters and it is generalised for both DM and MOND. 

The total model rotation speed is required from the total model 
potential in order to compare with the measured rotation curve. 
This can be found from 


VtotiRf 

R 




dR 


-(R). 


( 1 ) 


Here, R is the cylindrical radius in the disk mid-plane and is 
the total gravitational potential. Next, the total potential is required 
from the total mass distribution. This is computed via the Poisson 
equation, which in Newtonian dynamics is 


V^^n,tot{R, z) = ^nGptotiR, (2) 

where G is Newton’s gravitational constant, ptot{R,z) is the to¬ 
tal mass density from all sources (see the following and i]2.1.1t 
and ^n,tot{R, z) is the total Newtonian potential. In Newtonian 
dynamics, ^n,tot{R, z) is fully equivalent to $tot{R, z), but ptot 
consists of pbar and poM- In MOND, ptot is fully equivalent to 
Pbar, but a seco nd step is made to find ^tot from i&n.tot, which is 
dMilgromlbOlOh 


V^^tot{R,z)^V ■ [t/(|V<E-„ ,tot j/ao) Vthn.fotj , (3) 

where ag ~ 3.6( kms“^)^pc“^ is the acceleration threshold of 
MOND and 1 / is an interpolating function, chosen here to depend 
on its argument as 


i'j(v) 


l + (l + 4y-^ 
2 


1 / 2-1 1 /t 


(4) 


where 7 = 1 is the simple t^-function and 7 = 2 is the standard 
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^/-function l lFamaev & McGauehll2012l Eqs. 51 and 53), or likewise 
by 

..(.)= (5) 

2.1.1 Baryonic density 

The last question here is the distribution of ptot. The deriva¬ 
tion of the contribution from atomic and molecular gas is described 
in DMSvii. Both are assumed to have non-smooth, axisymmetric, 
radial surface densities. Both gas components are included in the 
modelling, hut they are considered fixed in mass and distribution. 
They are given nominal scale-heights of 200 pc, to which reason¬ 
able variations are inconsequential. 

There are typically two stellar components, a bulge and a disk. 
Below, the bulge and disk surface brightnesses are given to ex¬ 
pose the free parameters in their fitting. The bulge is assumed to 
be spherical and to follow a Sersic profile 

h{R) = h - 1]}, (6) 

where the effective surface brightness (7e), the effective projected 
radius (Re) and the Sersic index (n) can be fitted to the observed 
surface brightness distribution^ Simultaneously, the disk luminos¬ 
ity density can be expressed with a simple form 

where Ld is the total luminosity of the disk. Hr is the scale-length 
and hz is the scale-height. The vertical distribution is characterised 
by the ex ponential function, but the commonly used se cld(z/zo) 
function (Ivan der Kruit & Searlelll98l1 : lBottemalll993h would be 
equally appropriate. Here zq is the sech^ scale-height which corre¬ 
sponds to 2hz at large Using the secld vertical distribution does 
not change the conclusions (see El. 

In general, the surface density can be found by integrating 
along the line of sight. For a face-on galaxy, the luminosity den¬ 
sity of Eq|7]can be projected to give the surface brightness 

Here the density profiles are always assumed to be smooth. To find 
the mass density of the bulge and disk, the luminosity densities of 
the bulge and the disk must be multiplied by their respective M/L 
so that 

P* = Tbib + Tdjd, (9) 

where T b and T d are the stellar M/L values of the bulge and disk 
respectively. The total density of baryons is then phar = p* + pg, 
where pg is the atomic and molecular gas density. As stated previ¬ 
ously, in MOND, pbar is equivalent to the total mass density, ptot 
(EqO, because there is no DM in MOND galaxies. For Newtonian 
gravity ptot = Pbar + Pdm. 

For the baryonic mass models, there are 8 parameters: 
Tb,Ie,Re,n,Ld,'Td,hR,hz. Of those 9, the surface brightness 
parameters are either directly observed or unambiguously fitted to 
the surface brightness profile. This leaves only Tj,, Td, /iz. Since 

^ The constant k is fixed at 7.67. 


the DMS sample is chosen so that the total bulge luminosity to to¬ 
tal disk luminosity is low, Tt, is relatively insignificant and is never 
independently varied in the modelling performed here (Tt = T^). 
Therefore, T d and hz are the only two free parameters that are rel¬ 
evant to a theory like MOND. In principle, the inclination of the 
galaxy also has a small amount of freedom but it is strongly cur¬ 
tailed by the luminous TF relation. The aforementioned free param¬ 
eters are fitted for through a simultaneous comparison of the model 
vertical velocity dispersions and rotation curves with the observed 
ones, as is described in ( 13.31 

2.1.2 Inclination 

The derivation of the rotation curve of a moderate or 
high-inclination disk galaxy from the measu red 2D velocity 
field, permits the fittin g of tilted rings (see iBegemanI 1 19891 : 
Ivan der Hulst et aPl 19921) . These tilted rings allow us to model the 
variation in the inclination and position angle of numerous con¬ 
centric annuli at different galactocentric radii. The inclinations of 
the various rings, as a function of radius, can vary by 10° (e.g. 
Ide Blok et al.ll2008l) depending on the quality of the data, the regu¬ 
larity of the velocity field and characteristic inclination of the disk. 

The DMS galaxies have low inclinations (they are close to 
face-on), thus the inclination has a lot of leverage on the inferred 
rotation speed because of the shape of the sine function. It is pos¬ 
sible to derive accurate kinematic inclinations for nearly face-on 
disks using the tilted disk (a s opposed to tilted ring) technique of 
lAndersen & Bershadvl ( l2013l) . It is al so possible to inf er the inclina¬ 
tion using the luminous TF relation dVerheiienllTOOlh . This relates 
the absolute K-band magnitude of the galaxy, Mk, to a measure of 
the outer rotation speed, Vf such that 

Vf = 0.5 X 10(512-Mk)/11.3 (10) 

Anderson et al. (2015, in prep) have determined that the kinematic 
and TF inclinations for the DiskMass galaxies generally agree well, 
although there are some outliers. 

Relating the measured, inclined outer rotation velocity 
VobsSinif) with the expected outer velocity from the TF relation 
(Eg I lot allows the expected inclination to be deduced. This incli¬ 
nation is only that expected for the outer parts of the rotation curve, 
and thus the inner parts can vary somewhat due to a warp. 

In addi tion to the Luminous TF relation, there is a baryonic 
TF relation dMcGaugh et al.ll20o5 : lMcGaugl]|2005l) which relates 
the total baryonic mass of a galaxy to its outer, flat rotation speed, 
where Vf oc GMb ar • 

This relation is fundamental to MOND and the exponent e = 4 
and the constant of proportiona lity ap are predic tions which agrees 
well with the observed relation dMcGaughl2005l) . Thus, the MOND 
baryonic TF relation can be written in a similar form to EgllOlas 

I/'* 

^ = Mg + (Tbfb+Tdfd) X io(^k,0-Mr)/2.5^ 

Gao 

Here, fb = fd are the fractions of the total luminosity 

contributed by the bulge and disk respectively. The absolute mag¬ 
nitud e of the Sun in the K-band is Mk,q = 3.28 dBlanton et al.l 
l2003h and Mg is the gas mass. 

It was found that the inclinations from Eq[TT]are typically be¬ 
tween 5-15% larger than those found with EqjTO] depending on the 
M/L used (here the M/L was taken to be between 0.6 and 1): 
a smaller M/L implies a smaller corresponding rotation velocity, 
hence a larger inclination. 
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For a large enough sample there should, in principle, be no 
systematic deviation from either TF relation. When modelling ro¬ 
tation curves in general, it is not always clear when the rotation 
curve has reached the terminal velocity, so some margin of error 
must be granted. Since Eq[T0]has no dependence on M/L, and to 
make it easier to compare with the DMS results, the luminous TF 
relation inclinations (Ea llOt are used in this article. 


2.2 Stellar vertical velocity dispersions 

2.2.1 Choice of stellar velocity dispersion ellipsoid parameters 

In addition to the measured rotation curve, the DMS also mea¬ 
sured the line of sight velocity dispersion profile of the stars over 
the full projected area of the disk. This is then azimuthally aver¬ 
aged to give a ID line of sight velocity dispersion, which can be 
conv erted to a vertical v elocity dispersion (af) through the equa¬ 
tion iWestfall et al ]l^ 


2 


cos^ i 


1 + 


tan^ i 
2a2 


(1 + / 3 ') 


-1 


( 12 ) 


Here, the inclination of the galaxy to the line of sight is again, i, 
and a = — &. B = — provide information about the stellar ve- 
locity dispersion ellipsoid. Generally, a and B are expected to take 
on certain values from measurements in the Solar neighbour hood 
dSinnev & Merrifiei3l 19981 : Toerssen & Shapiro Griffir]]2oI^ . and 

B is presumed to take on specific values from the epicycle approxi¬ 
mation. However, beyond the Milky Way their variation is not em- 
piric ally well known (see discussion in sect 2.1 of DMSii , and 
also IWestfall et al.ll201ll: [WestfalllbOli iGentile et alj|2015h . By 
choosing galaxies that are nearly face-on, the DMS reduces their 
importance (cf. Eig|3- The statistical variation of these parame¬ 
ters is discussed in DMSii §2.1, and the DMS analysis establishes 
Q = 0.6 ± 0.15 and B = 0.7 ± 0.04. The mean values used by the 
DMS are chosen for the default values in this article. 


2.2.2 Model Stellar vertical velocity dispersions 

In order to compare with the observations, the model vertical 
velocity dispersions of the galaxies must be computed. The vertical 
velocity dispersion at a height z above the mid-plane, at a radial 
distance R from t he centre of the disk galaxy is found from (see 
iNipoti et'n]|2007l) 

P:t{R,z)az{R,z) = J p*{R,z) - — - dz , (13) 

and the equivalent of the observed vertical velocity dispersion at 
any radius, R, weighted by the local stellar surface density is given 
by 


/ OO 

pt{R,z)az{R,z)'^dz. (14) 

- OO 


These equations effectively reduce to 


•1 POO POO 

dRf = J J exp{-z'/hz) 


d^tot{R,z') , 


dz. 
(15) 

In the Newtonian gravity framework. Eg 1 151 depends mainly 


on 2 fitted parameters: Td and h^. As stated previously, inclina¬ 
tion could also be varied but only in a tight range around the TF 
relation values. Thus, combining simultaneous fits to the observed 
rotation curves and vertical velocity dispersions is a strong test of 
the MOND paradigm. 

It is worth noting that the infinite potential well of isolated 
galaxies in MOND is irrelevant here since the vertical gravitational 
field in EqlJ^is convolved with the exponentially declining stellar 
density, and thus the MOND gravity is only relevant where there 
are stars. 


2.2.3 The DiskMass Survey Method 

In the analysis of DMSvi, the measured rotation curve of each 
galaxy is used to fit the DM halo parameters and then the vertical 
velocity dispersion is used to directly give the mass surface density 
using T,dyn{R) = , where k is assumed to be 1.5 to describe 

an exponential vertical stellar distribution. From this surface den¬ 
sity, the gas disk was subtracted. This left the stellar disk surface 
density, E*(i?) = "EdyniR) — T,gas{R). This also includes an un¬ 
known component of DM. The stellar M/L as a function of radius 
was given by Td(i?) = An average of this M/L out to a 

given radius then defines the quoted disk M/L. This approach, al¬ 
though straight-forward, is only accurate when the derivative of the 
rotation curve is small. It is also not transferable to MOND because 
of the non-linearity of the theory and that the M/L in MOND af¬ 
fects both the rotation curve and vertical velocity dispersion. It is 
more secure to make the reverse calculation and go from observed 
surface brightness, sample a M/L to give surface density, then use 
EqUSjto give a model vertical velocity dispersion which can be 
compared with the observed vertical velocity dispersions. However, 
this requires calculation of the full three dimensional potential from 
a model galaxy - which is performed here. 


3 PRELIMINARY MODELLING 
3.1 Bulge-Disk decomposition 

The inclination corrected surface brightness profiles presented 
by DMSvi were analysed with a Markov Chain Monte Carlo 
(MCMC) approach to fit bulge plus disk surface brightness mod¬ 
els using Eqs[^and[8l The seeing was simultaneously accounted 
for with a Gaussian convolution of appropriate radius (see tables 3 
and 4 of DMSvii). All 5 surface brightness parameters (2 disk and 
3 bulge) were fitted, but the bulge surface brightness parameters 
were found to be degenerate with each other due to seeing effects 
and the lack of data points at low radii. Thus, the bulge parameters 
were fixed at the maximum likelihood values, which allowed better 
sampling of the more important disk parameters. In Figs ll Ill'S] the 
posterior probability distributions are presented for the two disk 
surface brightness parameters of each galaxy. For the disk scale- 
length, the range found by DMSvi (red curve) is also plotted. 

The best fits to the surface brightness profiles are given in 
Figs 1 1411^ and are generally good. Both linear radius and log ra¬ 
dius are plotted on the x-axis to expose the quality of the fits to the 
bulge and the outer disk. The fitted scale-lengths found here (see 
Figs ll llOt are generally consistent with those fitted by the DMS, 
although theirs have smaller error ranges. 

Once the bulge and disk surface brightness profiles were es¬ 
tablished, an Abel transform was used to de-project the bulge sur¬ 
face brightness, which was stored numerically. The disk has cylin- 
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Figure 1. Vertical gravity profile through three discrete cylindrical radii in the disk mid-plane R = 3.5, 8.5 and 18.5 kpc. (top panels) Comparison of the 
Newtonian theoretical vertical gravity profile against the computations from a Poisson solver for two different model galaxies. The galaxies are composed, 
as per Eq|7] of a radial exponential disk and exponential scale-height. The left hand panels use a galaxy with scale-length hji = 2.5 kpc, scale-height 
hz = 0.3 kpc and total mass Md = The right hand panels use a galaxy with scale-length hn = 4 kpc, scale-height hz = 0.4 kpc and total mass 

Md = 5 X The computations from the Poisson solver with Newtonian dynamics use the thick black line and the theoretical values found using 

Kuijken & Gilmore (1989) Eq 27 use the thin red line. The short blue vertical lines mark the scale-height of the model disk. The match is generally very good 
for 2 > O.lfcpc. (bottom panels) Comparison of the MOND (using 7 = 1; turquoise lines) and Newtonian (black lines) vertical gravity profiles as calculated 
by the Poisson solver. The galaxies used in both panels are the same as the ones used in the top panels. In these units, the MOND acceleration parameter is 
ao = 3.6 ( kms“^)^pc“^. Clearly the MOND boost to the gravity is relatively less significant for the higher surface density galaxy. 


drical symmetry, so it does not require de-projection. Following 
this, rejection sampling was used to generate N-particle represen¬ 
tations of the galaxies, with half of all particles representing the 
stellar disk and the other half split evenly between the stellar bulge 
and the atomic and molecular gas disks. The fitted scale-lengths are 
used for all the following analysis, but they are not used to generate 
the N-particle representations of the stellar disk. Instead, the “ob¬ 
served” surface brightness of the disk is sampled after subtracting 


the fitted bulge - since the fit to the surface brightness can be poor 
at large radii. 

In this analysis the static N-particle representations are only 
used to find the potential and gravitational field, there is no N- 
body evolution of the simulated galaxies. The masses recovered 
after generating the N-particle representations of each component 
of each galaxy generally agreed very well with the masses reported 
by the DMS. 
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3.2 The MOND Poisson solver 

The MOND Poisson solver described in I Angus et alj 1 I 2 OI 2 I 1 
is used to compute the radial and vertical gravitational fields of the 
N-particle galaxy models. The code solves the modified Poisson 
equation of the Quasi-linear version of MOND (QUMOND) given 
by Eq|3 The code uses a 3D grid and the cloud-in-cell technique to 
numerically discretise the 3D density of an N-particle distribution. 
It then employs finite differencing and multigrid methods to iterate 
from a test potential to the final potential which accurately reflects 
the density. 


3.2.1 Comparison of theoretical and numerical vertical gravity 

The Poisson solver is used to compute the gravitational field 
of the baryons in the radial direction, (R), as a function of 

radius, to find the model rotation speed (EqQ}. Simultaneously, the 
gravitational field of the baryons is solved for in the vertical di¬ 
rection, ^ gg ^ function of height above the disk, at sev¬ 

eral discrete radii: R = 0.5 to 9.5 kpc in steps of 1 kpc and then 
R — 11.5, 14.5, 18.5 and 25.5 kpc. In MOND, both gravitational 
fields noted above are equivalent to the gradient of the total poten¬ 
tial (i.e. ^bar = 4?tot) used in EQsfTlandlOland thus the rotation 
speed can be straight-forwardly calculated and Easll3landll4lcan 
be integrated to find at different radii, R. 

The most important, and difficult to compute, quantity is the 
vertical gravity profile at large radii. For an isolated, double expo¬ 
nential disk, like that introduced in Eq [T] the Newtonian vertical 
gravity profile at a gi ven radius can be ca l culate d numerically - as 
described in detail bv iKuiiken & Gilmorel (Il989h . 

In the top panel of Fig [T] the Newtonian vertical gravity pro¬ 
file is plotted for three discrete cylindrical radii R = 3.5, 8.5 and 
18.5 kpc. The left and right hand panels correspond to two differ¬ 
ent double exponential disk galaxies (as per EqO. The left hand 
panel has scale-length hu = 2.5 kpc, scale-height hz = 0.3 kpc 
and total mass Md = Mq. The right hand panel has scale- 
length /i_R = 4 kpc, scale-height hz — 0.4 kpc and total mass 
Md = 5 X 10^'^ Mq. The computation from the Poisson solver us¬ 
ing Newtonian dyna mics is the thick black line and the theoretical 
value using Eq 27 of iKuiiken & Gilmorel ( 1 19891) is the red line. The 
agreement is generally very good, except for small heights above 
the disk. This is due to a combination of the limit of spatial resolu¬ 
tion at large radii, owing to the centrally refining mesh, and limited 
particles at large radii. The average percentage error for these dis¬ 
crete radii (3.5, 8.5, 11.5, 14.5, 18.5 and 25.5 kpc) between the 
numerical calculation of Uz (using Eq[T5j and the theoretical val¬ 
ues are 0.05, 0.8, 0.7, 5.8, 5.2 and 4.4%. Given that only 6 of the 30 
galaxies have vertical velocity dispersion data points beyond 11 kpc 
and the errors on those data points are typically more than 10 %, this 
is more than adequate. 

In the bottom two panels of Fig[T]the Newtonian (black lines) 
and MOND (7 = 1; turquoise lines) vertical gravity profiles are 
compared, both of which were computed with the Poisson solver. 
The MOND profiles are strongly boosted relative to the Newto¬ 
nian profiles, however, the higher surface density galaxy (left hand 
panel) receives less of a boost than the other. 


3.3 Error budget 

The extra information that allows the DMS to close their set 
of equations is that, in general, disk scale-heights are observed to 


correlate with their scale-lengths. These scale-heights and scale- 
lengths have been fitted to edge-on galaxies and certain assump¬ 
tions are made in their modelling (like constant inclination, that 
the disks are well described by exponential or sech^ distributions). 
The DMS made an analysis of literature measurements and derived 
a simple relationship between hn and hz such that in units of kpc 
(see DMSii) 

hz ^0.2h°jf^^. (16) 

This relation has a Icr scatter of roughly 25% (DMSii). 

In the analysis presented here, the scale-height and disk M/L 
are fitted to the observed vertical velocity dispersions and rotation 
curves of the sample of 30 galaxies. After the distribution of fit¬ 
ted scale-height and scale-lengths is known, they will be compared 
for consistency with the direct observations of hu and hz for the 
sample of edge on galaxies compiled by the DMS. Therefore, a 
simultaneous fit must be made to the observed vertical velocity dis¬ 
persions and rotation curves. A critical concern is the contribution 
of each of the two datasets to the overall likelihood. Since the verti¬ 
cal velocity dispersions are the primary data set, the errors on each 
data point as taken as computed by the DMS. Given that there is the 
potential for some systematic errors from disk warping, the rotation 
curve should not be given too large a weighting, especially the in¬ 
ner parts. The decision was made to increase the error bars on the 
rotation curve data points to 10 kms“^. Lastly, the two separate 
reduced y^s from fitting the vertical velocity dispersion profile and 
the rotation curve are combined. 

Therefore, the likelihood is given by 
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where urc and nvvD are the number of relevant data points 
in the rotation curve and vertical velocity dispersion profile respec¬ 
tively. The prior then multiplies C to give the un-normalised poste¬ 
rior probability. 


4 PRIMARY RESULTS 

4.1 Parameter posterior probabilities 

MCMC sampling was used to find the posterior probability 
distribution for each of the free parameters: disk scale-height (hz), 
disk M/L (Td) and inclination (i). The scale-height was varied by 
considering the ratio of the fitted scale-height to the one derived 
from observations (Eq [HJ. A broad Gaussian prior of width 1.5 
was placed on this ratio. The prior placed on the disk M/L was 
centred on Td = 0.3 Mq/Lq and had a width 0.5 Mq/Lq. The 
ratio between the fitted inclination and the luminous TF relation in¬ 
clination (from EallOt received a fairly tight Gaussian prior of 0.15 
(or 15%), since there should be little deviation from the TF relation. 
The main reason for including inclination as a free parameter is just 
in case there is a sharp increase in posterior probability for a small 
change in inclination. The full expression for the likelihood is given 
by Eq[T3 This defines how the goodness of fit for the two data sets 
is combined: rotation curve and vertical velocity dispersion. 

The re-normalised posterior probability distributions of each 
of the three parameters (scale-height, inclination and M/L) is plot¬ 
ted in Figs ll7ll^ for each galaxy individually. There are two main 
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lines in each of the three columns: MOND with 7 = 1 and 2 (see 
Eq|4} black lines, solid and dashed respectively. The second col¬ 
umn (inclination) also has the curve of the prior (turquoise). The 
third column (M/L) has a green line which represents the Newto¬ 
nian gravity (with DM halo) MjL found by the DMS (DMSvi), 
shown only for reference. All the lines are re-normalised to have 
the same maximum posterior probability. The two MOND fits vary 
little in terms of goodness of fit. 

From the left hand column it is clear that the majority of 
the fits require substantially lower scale-heights than those derived 
from observations (Eq|16b. which in those panels are unity. There is 
a preference with most galaxies to have a higher inclination because 
that decreases the amplitude of the rotation curve. A lower rotation 
curve requires a lower M/L. A lower M/L allows a larger scale- 
height - which allows better agreement with the observations of 
edge-on galaxies (Eollbb. However, the fairly tight prior on the in¬ 
clination prevents it from changing substantially. Usually it changes 
by less than 10 %. 

In Table [T] various Icr confidence ranges for the fitted scale- 
heights and stellar mass-to-light ratios are presented using MOND 
with varying constraints (from interpolating functions and which 
parameters are left free). The preferred scale-heights, those from 
Eq[T^and used by the DMS, are given for each galaxy (column 2). 
For reference (column 3) the M/L found by the DMS when using 
Newtonian gravity and DM halos is shown, followed by the best 
fit MOND scale-heights (as a ratio between the fitted scale-height 
and the one derived from observations of edge-on galaxies - Eq[T 6 ) 
and the M/L for three different interpolating functions (the two 
used above and one other). These fits are for the scenario where the 
scale-height and M/L are barely constrained, but the inclination is 
tightly constrained by its prior. 


4.1.1 Comparison of the MOND fits with the DMS data 

In Figs 123125 ] the best fits to the vertical velocity dispersion 
(first and third columns) and rotation curve (second and fourth 
columns) are plotted in the MOND case with 7 = 1 . The two inter¬ 
polating functions are not plotted together because different best fit 
inclinations mean the data points vary in each case. There is little 
difference between the two MOND cases in terms of quality of fit. 
Most fits to the vertical velocity dispersion are good, but the qual¬ 
ity of fits to the rotation curves range from poor to good. Rotation 
curves with good fits include the galaxies UGC 448, 1081, 1087, 
1908, 3091, 4036, 4368, 4380, 7244, 6918, 8196, 9965, 11318, 
12391. Quite good fits are found for UGC 463, 1529, 1635, 1862, 
3701, 4256; whereas the fits to UGC 3140, 3997,4107, 4458, 4555, 
4622, 6903, 7917, 9177, 9837 are poor. 

MOND fits to the rotation curves of n early face-on galaxies 
are ve ry sensitive to the inclination (see e.g. Ide Blok & McGaurf] 
Il998h . In the outer parts, they depend on the fourth power of 
l/sin{i). Warping of the disk might be invoked to improve cer¬ 
tain rotation curve fits, however, in some cases the DMS galaxies 
have already been corrected for warps. 

To compare with models that use the scale-height from Eall6l 
the vertical velocity dispersion (red line) is over-plotted for each 
galaxy using that larger scale-height, but the M/L and inclination 
of the best fit are kept. In many cases, the red line has a far greater 
amplitude than the data points and best fit (black line). There are, 
however, a few cases where the red and black lines are close to¬ 
gether, such as UGC 1862, UGC 4368 and UGC 4458 - although 
the latter only has a single relevant data point. 


4.2 Comparison of the fitted hz vs ha with observations 

From fitting the observed surface brightness profiles of the 
DMS galaxies, a value for the radial scale-length was estimated 
for each galaxy disk. As can be seen in Figs II 1113] (second and 
fourth columns) the radial scale-length fitted by DMSvi is gener¬ 
ally slightly larger than the best fit found here, but the majority are 
consistent within the errors. This fact makes the analysis presented 
here slightly more favourable to MOND. The difference in fitted 
scale-lengths is in part due to the fact that in DMSvi a finite central 
region was used to fit the scale-length, whereas the full extent of 
the galaxy is used here. Thus, for any galaxies whose disk surface 
brightness is not well fit by a single exponential there is a differ¬ 
ence. For example, UGC 4036 has a shallow inner surface bright¬ 
ness profile and then a steep outer decline, thus the scale-length 
found here naturally has a smaller scale-length. Recall that the N- 
particle realisations ( i|3.1b use the observed surface brightness, not 
the fitted exponential disk. 

Similarly, by fitting the vertical velocity dispersion profile and 
the rotation curve of each DMS galaxy, a confidence range for the 
scale-height of each galaxy has been derived. In principle, this com¬ 
bination of hz and ha parameters can be compared to the measure¬ 
ments of scale-lengths and scale-heights of a sample of edge-on 
galaxies. In Figl^the hz vs ha diagram for the MOND fits (7 = 1, 
ao ^ 3.6( kms“^)^pc“^) is plotted along with the observations 
of these parameters for edge-on galaxies made by K02. 

In order to visualise this more plainly, contours in the hz vs 
ha plane were added to Fig|^ Here, the re-normalised posterior 
probability density of all individual data points are co-added on a 
fine grid, with each galaxy receiving equal weighting. A point with 
larger error bars will spread its probability density over a larger 
range. The blue contours represent the co-added MOND probabil¬ 
ity densities and the black triangles give the location of each in¬ 
dividual galaxy point. The error bars are not included to minimise 
clutter. The measurements by K02 have red shaded contours with 
green circles. The MOND contours are significantly inclined rela¬ 
tive to the K02 contours. The equivalent plot for the 7 = 2 interpo¬ 
lating function is almost identical. 


4.3 Two dimensional Kolmogorov-Smirnov test 

To demonstrate statistically the the large offset in the hz vs ha 
plane between the K02 data points in Fig [2] and the fitted MOND 
points, a 2D Kolmogorov-Smimov (KS) test (see jPress et al ■I 19921) 
chapter 14.7) was employed. In this test, the first K02 data point 
was taken and made the origin of four quadrants. Then the fraction 
of K02 points in each quadrant and the fraction of MOND points 
in the same quadrants. The largest difference between the two frac¬ 
tions in a single quadrant is stored and then the rest of the K02 
points are cycled through, making each one the origin in turn. The 
largest of the largest differences is stored and then the procedure 
is repeated centring on the MOND points. The largest difference 
in fractions is ordinarily the same regardless of whether the K02 
or MOND points are central. This number is the statistic of merit 
for the 2D KS test and provides a significance level that the two 
samples originate from the same parent population. 

Since the data points have fairly large error bars, if the mean 
of a particular data point lies in the first quadrant, it is not ideal to 
add the whole point to that quadrant. Instead, the fraction of the 
point inside each quadrant determined by the point’s error bars is 
added. 

According to the 2D KS test, the significance levels that the 












8 G.W. Angus, G. Gentile, R. A. Swaters, B. Famaey, A. Diaferio, S. S. McGaugh, K. J. van der Heyden 


1.0 


^ 0-8 
o 

Qh 

0.6 


X, 

QjO 

X 0.4 

0 ) 

"cfl 

o 

m 

0.2 


0.0 

0 2 4 6 8 10 

Scale length, hj; [kpc] 

Figure 2. Contours of the posterior probability distribution of galaxies in the vs hji plane. The blue contours and the black triangles correspond to the 
scale-lengths fitted to the surface brightness profiles of the DMS galaxies and the scale-heights found from the MOND (7 = 1, oq = 3.6( kms“^)^pc“^) 
fits to their vertical velocity dispersions. The measurements of scale-lengths and scale heights from observations of edge-on galaxies by K02 have red shaded 
contours with green circles. The re-normalised contour levels are 1.0, 0.91, 0.82, 0.73, 0.64, 0.55, 0.46, 0.37, 0.28, 0.19, 0.1, 0.0001. The two sets of contours 
are highly inclined relative to each other. 



MOND fits with 7 = 1 and 2 come from the same parent distri¬ 
bution as the K02 sample are 4.4 x 10“® and 1.4 x 10“"^ respec¬ 
tively. Using the interpolating function <5 = 4 the significance level 
is 1.0 X 10“^. Comparing with DM models, it is easy to obtain 
significance levels of order unity, so the method seems robust. 

As a separate test, a straight line was fitted to each data set of 
hz vs /i_R with a fixed intercept. The resulting gradients differ at 
the 2a level, regardless of interpolating function. Using a variable 
intercept only improves the agreement slightly. 


4.4 Superposition of M jL distributions 

In the approach presented here to modelling the DMS data 
the MjL is fitted for as well as the scale-height. At near infrared 
wavelengths a smaller variation of the MjL is ex pected from one 
galaxy to the next th an at optical wavelengths teell & de Jond 
I 2 OOII : iBell et al.]|2007h . To find the posterior probability distribu¬ 
tion of MjL across the full sample of galaxies, the summation of 
the re-normalised posterior probability density for the MjL of all 
30 galaxies is plotted in Fig[^ This is done for three MOND inter¬ 
polating functions and the DMS computation of MjL in the DM 
scenario. This means the same weighting is given to each galaxy’s 
MjL posterior probability density in the summation. Recall that in 
the fitting process a Gaussian prior is used with width O.SMq/Lq 
centred on M/L = O.SMq/Lq. 

The 7=1 MOND interpolating function (solid black line) 
has an average best fit M/L ~ 0.55 ± 0.15, whereas the 7 = 2 
MOND interpolating function (dashed black line) requires signifi¬ 
cantly larger values, M/L ~ 0.7±0.2. Bearinmindthatthe 7 = 2 
produces less of a boost to the gravity than the 7 = 1 interpolating 
function at intermediate gravities around ~ oq. The <5 = 4 interpo¬ 
lating function (which uses Eq[^ has a similar M/L distribution 
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Figure 3. Co-addition of all 30 re-nomialised posterior probability profiles 
for the fitted M/L parameter. MOND fits for the interpolating functions 
using Eq|4]with 7 = 1 and 2 are represented by black solid and dashed 
curves respectively and 5 = 4 (from Eq|^ has a dotted curve. The co-added 
M/L found by the DiskMass Survey in DMSvi has a green curve. 

as 7 = 1. There is no impact on the derived M/L values from the 
weak prior. The green line, given only for reference, is the M/L 
derived from the DMS analysis with Newtonian gravity and DM. 

4.5 Why does MOND need such thin disks? 

The results presented here suggest that MOND requires disks 
that have roughly half the vertical scale-height as those inferred 
from observations of edge on galaxies (DMS II or Eg I16t. There 
are two key differences between the MOND vertical velocity dis¬ 
persions and those using DM. The first effect is related to the re¬ 
quired M/L of the galaxies. Comparing a disk with the same sur¬ 
face brightness in MOND and Newtonian gravity, the Newtonian 
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disk is surrounded by a dark matter halo. This dark matter halo can 
account for the radial gravitation field that produces the observed 
rotation speed. This means the M/L of th e Newtonian disk is onl y 
weakly constrained by the rotation curve jvan Albada et al.lll985h . 
In MOND, the M/L is essentially fixed by the rotation speed be¬ 
cause it is the only parameter that can be varied to allow a match 
between the observed and expected rotation speeds. Based on the 
results presented here, the required M/L of the disk is a factor of 
between 1.5 and 2.5 higher in MOND (cf. Fig[3, depending on the 
interpolating function used. Since the M/L is typically much larger 
in MOND, the resulting vertical gravitational field is larger than the 
dark matter equivalent. Moreover, even if the disks in MOND and 
Newtonian gravity had the same M/L, the MOND gravitational 
field is still boosted relative to the Newtonian one as seen in Fig[T] 
This follows because the nearly spherical dark matter halo typi¬ 
cally has a weak influence on the vertical gravitational field within 
the disk. 

The influence of MOND on the vertical gravitational field, rel¬ 
ative to Newton, for the same disk is of course greater when the sur¬ 
face density is low. This is demonstrated by comparing the vertical 
gravitational field of a high surface brightness galaxy (Fig[T]bottom 
left panel) with a lower surface brightness one (Fig[T]bottom right 
panel). 

In relation to the equations of jJH it can be seen that increas¬ 
ing T d or hz increases the amplitude of the model However, 
in MOND only increasing Td increases the model rotation curve 
at intermediate to large radii. Thus, once is fixed by the rota¬ 
tion curve, only hz can be varied to fit the observed vertical ve¬ 
locity dispersion. Since the observed vertical velocity dispersions 
are typically much lower than predicted by MOND, hz must be de¬ 
creased to fit them. It is for these reasons, the fitted scale-heights in 
MOND are significantly lower than those expected from observa¬ 
tions (EqlUll. 

4.6 Central dynamical surface density versus central surface 
luminosity density 

Is waters et al.l ( |2014|) used the DMS data to present a correla¬ 
tion between the extrapolated disk central surface luminosity den¬ 
sity, from the photometry, and the extrapolated disk dynamical cen¬ 
tral surface mass density (see 112.2.3b . The majority of galaxies are 
consistent with the relation T,dyniO) = O.SMq/L©/,( 0), but at 
the low surface luminosity density end a clump of 6 galaxies lie 
above the relation. In the Newtonian context, this can be attributed 
to greater rela tive DM dominance for low surface luminosity den¬ 
sity galaxies dSwaters et al.ll20l3) . but this is also reminiscent of 
the MOND transition in rotation curves of spiral galaxies, where 
the effect becomes pronounced at low surface densities. Thus, a 
prediction for dynamical surface density versus central surface lu¬ 
minosity density in MOND would be valuable. 

To this end, a representative galaxy from the DMS was se¬ 
lected - UGC 4036. The stellar bulge and all gas was removed 
from the N-particle realisation, thus it has the same stellar disk sur¬ 
face brightness profile as UGC 4036, and the same vertical profile. 
For this model, the central vertical velocity dispersion, (0) from 
Eg 1151 can be found for a chosen M/L. From here, the surface 
brightness of the stellar disk can be scaled to compute az (0) at a 
broad range of central surface luminosity densities which allows a 
comparison with the data from is waters et al] l l2014l) . which is done 
inFig|4] 

The data points come from ISwaters et alj ( 1201 4h Fig 1 (left 
hand panel) and both coordinates are extrapolated from exponential 
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Figure 4. Disk central dynamical surface density versus disk central sur¬ 
face luminosity density relation. The data points con'espond to each galaxy 
in the DiskMass Survey and are reproduced from Swaters et al. (2014). The 
red line illusttates T,dyn{0) = 0.3Mq/LqI,{0). The other lines show 
MOND predictions: the dashed lines use scale-heights that are twice as thin 
as those derived from fits to observations of edge-on galaxies (Ea ll6l . and 
the solid black line uses the scale-heights derived from fits to observations 
of edge-on galaxies. The black lines all use M/L = 0.6. The upper blue 
dashed line uses M/L = 0.9 and the lower one uses M/L = 0.3. The 
appai'ent flattening of the central velocity dispersions at low surface bright¬ 
nesses is consistent with MOND. 

fits to the observed projected disk surface luminosity density and 
observed vertical velocity dispersion. The red line is the relation 
'Ldynfi) = 0.3Mq/Lq7*(0). The black solid line is the MOND 
prediction with M/L = 0.6 and the scale-height expected from 
observations of edge-on galaxies (DMSii; EqlJ^in this paper). The 
black dashed line has the same M/L, but used disks that are twice 
as thin, which MOND has been shown here to require in order to 
mesh with the DMS data. The two blue dashed lines also use disks 
twice as thin as observations of edge-on galaxies, with the upper 
line using M/L = 0.9 and the lower line using M/L = 0.3. 
Thus, MOND provides a natural explanation for the deviation of 
low surface luminosity density galaxies from the central dynamical 
surface density versus central surface luminosity density relation. 

5 SECONDARY RESULTS: VARYING ADDITIONAL 
PARAMETERS 

In this section the impact of other variables on the results from 
ll4lis considered. 

5.1 Variation of the Stellar Velocity dispersion Ellipsoid 

As discussed in 112.21 the stellar velocity dispersion ellipsoid 
(SVE) parameters are fixed to those used by the DMS i.e. a = 0.6 
and j3 = 0.7. To demonstrate the degree to which a reasonable 
variation in these parameters might affect the results, in Fig[5]the 
vertical velocity dispersion ratio between three combinations of pa¬ 
rameters and the default values are plotted against inclination. 

Using different values for a in EqlJ^can make a sizeable dif¬ 
ference to the derived Gz, especially for large inclinations. The dot¬ 
ted line of Fig|5]shows that increasing a to 0.9, which is 2a above 
the default value used by DMSvi, would increase the derived Gz 
by around 10% for standard DMS inclinations of 25°. On the other 
hand, decreasing a to 0.48 and using /3 = 1.04 (dot-dashed line 
in Fig |5]l, which was the value derived by IWestfall et alj ||201 ih 
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Figure 5. Ratio between the Uz found from Ea ll2l using certain a and /9 
parameter combinations, and the default values used by the DMS (which are 
a = 0.6 and /3 = 0.7 as per DMSvi) against inclination. The combination 
of a and /3 used for each line is made clear in the hgure. Varying a can 
have a significant impact on the derived velocity dispersion for increasingly 
inclined galaxies. 


for UGC 463, would decrease the derived Gz by around 15% for 
i = 25°. 

To test the impact of these SVE parameters on the fits to 
the vertical velocity dispersions and rotation curves, MCMC mod¬ 
elling was used to simultaneously fit M/L, a and /3, with incli¬ 
nation fixed by the TF relation and hu/hz set by observations of 
edge-on galaxies (Eall6t. Flere, only flat priors are used such that 
0 < a < 1 and 0 < /3 < 2.5. The M/L has a Gaussian prior of 
width 0.5 centred on 0.3. The most relevant parameter is a, due to 
its impact on Eall2l 

In the top panel of Fig the co-addition (for all 30 galaxies) 
of the re-normalised posterior probability is plotted as a function 
of a. Each of the 30 galaxies are given equal weighting, regardless 
of their individual posterior probability. This shows the impact of 
increasing a from the nominal value of 0.6 used here and by the 
DMS. The posterior probability increases roughly a factor of two 
and a half from a = 0.6 to 1. 

This increase does not necessarily mean that the fit qualities 
are good. To demonstrate this, it is worth comparing the fits where 
a is a free parameter (but scale-height is fixed) to the fits where the 
scale-height is free. These fits where scale-height is free represent, 
for the majority of galaxies, a good quality fit to the vertical velocity 
dispersion. For the remainder of this section, this fit is referred to as 
the benchmark fit. a does not influence the fit to the rotation curve. 

In the bottom panel of Fig |6] the y-axis is the ratio between 
the posterior probability of the fit with fixed scale-height (Eg 116b 
and inclination (luminous TF relation), a = 1 and best fit /3 and 
M/L to the benchmark fit. This is plotted against the inclination of 
the galaxy derived from the luminous TF relation to show that in¬ 
creasing a may help with some (but not all) of the high inclination 
galaxies, but the majority of the sample would have extremely poor 
fits to the DMS data even with a — 1. 
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Figure 6. (top panel) Superposition (for all 30 galaxies) of the re¬ 
normalised posterior probability density as a function of a. (bottom panel) 
Ratio between the posterior probability of the MOND fit to the DMS data 
with fixed scale-height and inclination, but a = 1 (along with best fit /3 
and M/L), to the benchmark fit. This is plotted against the inclination of 
the galaxy derived from the luminous TF relation. The MOND parameters 
7 = 1 and ag = 3.6 (kms“^)^pc“^ ai'e used. 


5.2 Fixed scale-height, unconstrained inclination 

A relevant question is how far the fitted inclination would have 
to deviate from the luminous TF relation before a good fit to the 
DMS data could be achieved, whilst remaining consistent with the 
hz values derived from fits to the photometry of edge-on galax¬ 
ies. To do this, the scale-height was fixed to the aforementioned 
value (Eg 1 16b for each galaxy and fitted only for the inclination and 
M/L. The parameters 7 = 2 and ag = 3.6 (kms“^)^pc“^ were 
used. For the sake of brevity, in Fig |7] only the co-added, posterior 
probability densities for the two parameters are plotted. 

The top panel shows the M/L, which is shifted to values 
roughly three times smaller (cf. dashed line of Fig[3 than the orig¬ 
inal models with scale-height free. In the last column of Table [T] 
the M/L range is given for each galaxy in this scenario individ¬ 
ually. The M/L values are much smaller than the other models 
where inclination is constrained. In the bottom panel of Fig |7] the 
inclination, which is given as a ratio of the fitted value (if a) to the 
luminous TF relation value (irp), is plotted. With an unchanged 
inclination (ifa = ipp), the quality of the fits to the majority of 
the DMS galaxies are poor, but by increasing the fitted inclination 
to ifit/ipp = 1.4 the fits are of a high quality. This is partly owed 
to the resulting rotation speeds generally having lower amplitudes. 
Such significant deviations from VOTs luminous TF relation seems 
unlikely. 

To demonstrate this last point, in Fig the luminous TF re¬ 
lation from VOl (Eg 1 10b is plotted, which is used for the DMS 
and this study. The solid black line is the best fit TF relation and 
the black data points are the values of the DMS galaxies. The blue, 
green and red points are the positions of those galaxies if their incli¬ 
nations were increased by a factor 1.1, 1.2 or 1.4 respectively. The 
dashed and dotted lines are the la and 2a scatter in the observed 
TF relation. Clearly even changing the inclinations by a factor of 
1.1 would be in stark conflict with the TF relation. 
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Figure 7. Impact of varying inclination on the goodness of fit and the MjL. 
These models use fixed scale-heights lEg llht . but the inclination is left com¬ 
pletely unconstrained, (top panel) Co-added posterior probability density as 
a function of MjL, for all 30 galaxies, required for high quality fits to the 
DMS data when the inclination is left free, (bottom panel) Co-added pos¬ 
terior probability density against the ratio of the fitted inclination to the 
inclination derived from the luminous TF relation (VOl). The inclinations 
required for good fits to the DMS data are far larger than the inclinations 
derived from the luminous TF relation. 


5.3 Variable MOND acceleration parameter 

The DMS data seem to be at odds with MOND when using 
the default acceleration parameter ao = 3.6 (kms~^)^pc~^. To 
test this parameter’s influence, models were tested where the scale- 
height and inclination are fixed, but the MOND acceleration pa¬ 
rameter is variable - along with the MjL. Although varying the 
MOND acceleration parameter from galaxy to galaxy is not an ac¬ 
ceptable solution, this could identify a unique value for oo that is 
more consistent with the DMS data. In the top panel of Fig [^an¬ 
other co-added, posterior probability density is plotted, this time 
against the MOND acceleration parameter. This is done for three 
different interpolating functions (7 = 1, 4 and 8 - black, red and 
green lines respectively). The parameters 7 = 4 and 8 are trialled to 
see if the sharper transition from the MOND regime to the Newto¬ 
nian regime would affect the ability to fit the DMS data. It appears 
that little is gained overall by varying ao and the same is true on a 
case by case basis. This can be demonstrated by comparing the fit 
with ao free (and scale-height fixed) to the benchmark fit (similar 
to what was done in 1 )5. II and Fig [^bottom panel). In the bottom 
panel of Fig |3 for each galaxy individually, the ratio is plotted be¬ 
tween the posterior probability of the best fit using a variable ag 
(with fixed scale-height) and the benchmark fit. The vast majority 
of the galaxies would have poor fits with a variable ao and fixed 
scale-height. 

Varying ag suffers from the same drawbacks as varying the 
MjL. Increasing ao can allow a fit to the rotation curve with a 
lower M/L, but this simultaneously increases the vertical gravita¬ 
tional field, which in turn requires the scale-height to be reduced. 

In the second (with 7 = 8 ) and third last (with 7 = 4) 
columns of Table [T] the M/L range is displayed for each galaxy 
individually when the MOND acceleration parameter is variable. 


Figure 8. Luminous TF relation from VOl (’Fg llOi showing K-band mag¬ 
nitude against outer rotation velocity. The solid black line is the best fit 
luminous TF relation and the black data points represent the DMS galaxies. 
The blue, green and red points are the positions those galaxies would take if 
their inclinations were increased by a factor 1.1, 1.2 or 1.4 respectively. The 
dashed and dotted lines are the Itr and Icr scatter in the observed luminous 
TF relation. 


5.4 Sech^ vertical distributions 


Regarding whether the poor results were caused by us¬ 
ing an exponential function to describe the vertical scale-height, 
the DMS data was modell ed using a sech^{z/zo) function 
jvan der Kruit & SearlelfTosTI) to describe the vertical distribution 
of the stellar disk in the N-particle models. The scale-height used 
in the sech^ function is zg. In Fig[T0]the posterior probability den¬ 
sities are plotted for the ratios of the fitted ha to zg for the 7=1 
and 2 interpolating fun ctions . These are compared with da ta from 
iMosenko v et alj feoid) an d Bizvaev & Mitrono^ j2002l) . From 
Fig 5d of IMosenkov et al] ( | 201 (]|) the indicative range of observed 
hn/zg is between 12 and 22. A solid vertical red line is placed 
in Fig [To] to mark the ha/zg value where their histogram begins 
to decline and a dashed vertical red line to mark where their data 

effectively terminates. _ _ 

The 153 galaxies of iBizvaev & Mitronoval 1 I2OO2I) are also 
used to produce a probability density (green curve in Fig m 
to compare with the posterior probability density found here. 
The hit/zg distributions (solid and dashed black lin e s in F igllOt 
fitted here peak be y ond where IMosenkov et alj (1201 Ol) and 
IBizvaev & Mitronoval ( I2OO2I) have any evidence of such values. 
The distribution found here also extends to much larger values. 
To demonstrate the mismatch statistically, a 2D KS test was per¬ 
formed to_com£M£flie_/ifl_andzo_^stributions found here with the 
data of iBizvaev & Mitronoval ( I2OO2I) . The confidence levels found 
for the 7=1 and 2 interpolating functions were 6.0 x 10 ~^ and 
2.8 X 10“®. Evidently, the sech^ vertical distribution does little to 
aid the agreement between MOND and the DMS data. 


6 CONCLUSIONS AND DISCUSSION 

By assuming the scale-heights of face-on stellar disks are sim¬ 
ilar to the measured scale-heights of edge-on galaxies (DMSii), the 
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Figure 9. (top panel) Co-added re-normalised posterior probability density 
for all 30 galaxies, this time against the MOND acceleration parameter, 
ao- Three interpolating functions were used: 7 = 1, 4 and 8 (black, red 
and green lines respectively). The default value for the MOND acceleration 
parameter is ao = 3.6 (km s“^)^pc“^. (bottom panel) Ratio between the 
posterior probability of the best fit using a variable ao and the benchmark 
fit. The three interpolating functions are represented by the same colours as 
the lines in the top panel. Varying ao does not allow good fits to the DMS 
data. 


DMS found that stellar disks are sub-maximal, with K-band mass- 
to-light ratios of the order of M/L ~ 0.3 (DMSvi). 

In this paper the measured vertical velocity dispersions and ro¬ 
tation curves from DMSvi were analysed in the context of MOND. 
The problem is that the velocity dispersions measured by the 
DiskMass Survey can only allow for thick enough disks if the ver¬ 
tical gravitational field is weak. In Newtonian gravity this can be 
accommodated by lowering the M/L. In MOND, the M/L must 
be larger to account for the rotation curve in the absence of DM. 
This means the vertical gravitational field is strong. Moreover, it is 
amplified by MOND relative to the vertical gravitational field of a 
disk with the same M/L in Newtonian gravity. 

Thus, the most straight-forward way to reduce the model ver¬ 
tical velocity dispersions in MOND is to decrease the stellar scale- 
height. If these were reduced from the values derived from obser¬ 
vations of edge-on galaxies by roughly a factor of 2 then the DMS 
vertical velocity dispersions and rotation curves would be compat¬ 
ible with MOND. Regardless of stellar scale-length, these disks 
would have scale-heights between 200 and 400 pc. According to 
a two-dimensional K-S tests, such thin disks are strongly at odds 
with observations. 

All other key parameters involved in analysing the DMS data 
in MOND were modelled, such as the stellar velocity dispersion 
ellipsoid parameters, inclination, the choice of stellar vertical dis¬ 
tribution profile, the MOND acceleration parameter and interpolat¬ 
ing function. No clear way was found to reconcile MOND with the 
data. 

If the derived rotation velocities were lowered, by imposing 
30-40% higher inclinations, it became possible to simultaneously 
fit the low vertical velocity dispersions and low rotation speeds with 
low values for the stellar M/L. By design, since the outer rotation 
speeds are fitted, this is in accordance with the MONDian baryonic 
TF relation (with logarithmic slope fixed to 4, Eq II lb . Despite this. 



Mass-to-light ratio, W/L 



Figure 10. (top panel) Superposition of the re-normalised posterior prob¬ 
ability densities of fitted mass-to-light ratios, using a sech^ function to 
describe the vertical distribution of the stellar disk, for the 30 galaxies in 
the DiskMass sample. The black solid curve corresponds to the 7 = 1 
interpolating function and the black dashed curve to 7 = 2. The distri¬ 
bution found by the DMS, with an exponential vertical distribution, uses a 
turquoise curve and is just shown for reference, (bottom panel) Superposi¬ 
tion of the re-normalised posterior probability density of fitted scale-length 
to vertical scale-height, using a sech^ distribution, for the 30 galaxies in 
the DiskMass sample. The black solid curve con'esponds to the 7 = 1 in¬ 
terpolating function and the black dashed curve to 7 = 2. The green curve 
is the probability density of the sample of 153 galaxies from Bizyaev and 
Mitronova (2002). The vertical solid red line represents the value of hij/zo 
at which the sample of Mosenkov et al. (2010) begins to decline from a 
previously flat value. The dashed red line represents where the probability 
density of Mosenkov et al. (2010) drops to zero. Using a sech^ distribution 
does not improve the ability of MOND to fit the DMS without decreased 
scale-heights. 


these circular velocities are inconsistent with the luminous TF re¬ 
lation which is quite well established for high-inclination galaxies. 
To this end, it would be interesting to make a fully self-consistent 
analysis that models the effect of the varying inclination on the two 
dimensional velocity fields and the surface brightness profiles. 

A prediction was made for the scaling of the central dynam¬ 
ical surface density versus central surface luminos ity density in 
MON D and compared with the DMS galaxies (see lSwaters et all 
l2014h . MOND correctly predicts that the low surface luminos¬ 
ity density galaxies deviate from the simple relation T,dyn{0) = 

0.3Mq/LqL{0). 

In ad dition to viewing MOND as a modification of the law 
of gravity, iMilgromI ( 1201 ll ) has maintained that a modification of 
inertia might be better positioned to explain all available galactic 
dynamics data. In such a theory, the orbital trajectories of particles 
would have different inertias meaning the interpolating function in 
the radial direction would differ from the vertical directions. This 
is an intriguing possibility since the problem with MOND outlined 
here is that for the radial force defined by a galaxy, the correspond¬ 
ing vertical force is too large. It is interesting therefore that the 
vertical velocity dispersions of disk galaxies must increase to agree 
with MOND, but the obser ved radial velocity dispersions of some 
dwarf spheroidal galaxies ( lAngus|[2008l : I Angus et al.ll2014h must 
decrease. This anecdotal evidence should be followed up with more 
rigorous study. 

Something worth further investigation is the impact 
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from super-thin disks, much t hinner than the observe d scale- 
heights represented here i Schechtman-Rook et iD I2012I: 
ISchechtman-Rook & Bershadvl I20l4l . on the measurement of 
vertical velocity dispersions by the DMS survey. Specifically, 
how prominent are super-thin disks in all disk galaxies, and 
are the vertical velocity dispersions that the DMS measures 
meaningfully influenced by them? Th e super-thin disk found by 
ISchechtman-Rook & Bershadvl 12014) was truncated at around 
3 kpc, suggesting the influence is minimal. As discussed in DMSi, 
it is expected that stars in or near the mid-plane diffuse over time. 
The youngest stars may contribute to the light because of the OB 
stars, but they contribute little to the absorption lines, because 
those OB stars have weak spectral features (except for H and He). 

It is possible to make even stronger claims about the status 
of MOND-like theories, as well as conventional dynamics, with 
certain desiderata. From a theoretical point of view, this would be 
a detailed study of the impact of the neglected cross-term (tilt of 
the velocity ellipsoid) in Eq[T5] From an observational point of 
view, a larger sample of edge-on galaxies (like K02) with near- 
infrared photometry, to more precisely confirm the correlation be¬ 
tween scale-length and scale-height. It would also be beneficial to 
have an increased sample of nearly face-on galaxies, with measured 
stellar velocity dispersions and rotation curves and with greater 
sensitivity to be able to measure the stellar velocity dispersions to 
larger radii. 
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Table 1. This table presents various scale-heights and mass-to-light ratios with Icr errors for each of the 30 galaxies in the DMS sample. Column (1) gives the galaxy name, (2) the scale-height used by the DMS, 
(3) has the M/L fitted by the DMS using Newtonian gravity and a DM halo. The remaining columns all pertain to MOND fits. Columns (4)&(5), (6)&{7) and (8)&(9) give the best fit ratio of fitted scale-height to 
the one used by the DMS and the fitted M/L for three distinct interpolating functions respectively. These models come from analyses where the M/L and scale-height are almost unconstrained, and the inclination 
has a limited amount of freedom. The next two columns (10) and (11) are for models where only the MOND acceleration parameter and M/L were free to vary (inclination and scale-height were fixed). These two 
columns use different interpolating functions and show the M/L of the best fit. The last column is the best fit M/L found when the inclination is unconstrained, but scale-height is fixed. 
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Figure 11. Re-normalised posterior probability profiles for the two disk surface brightness parameters {Eq[^ fitted to the DMS galaxy surface brightness 
profiles (DMSvii). In the first and third columns the disk central surface brightness is plotted and the disk scale-length is plotted in the second and fourth 
columns. For the disk scale-length, the values found by DMSvi (red curves) are over-plotted for comparison. 
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Figure 12. As per Fig[TT] 
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Figure 13. As per Fig[TT] 
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Figure 14. Obsei'ved surface brightness profiles {black error bai's around small circle) along the disk major axis for each galaxy in the DMS. Each row shows 
two different galaxies. In columns 1 and 3 the fits in log-log are plotted and in columns 2 and 4 they are plotted as log-linear, to expose the quality of the fits at 
different scales. The green dashed lines are the intrinsic bulge surface brightnesses and the solid green lines ai'e the seeing affected versions. The same is true 
for the red dashed and solid lines, except these are for the exponential disks. The black solid lines are the seeing affected total combined surface brightnesses. 
The vertical turquoise lines are the bulge radii as given by DMSvi and are unchanged in this analysis. 
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Figure 15. As per Fig[T5] 
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Figure 16. As per Fig[T3] 
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Figure 17. Re-normalised posterior probability for each parameter from the fits to the vertical velocity dispersions and rotation curves. Each row displays a different galaxy. The first column shows the ratio of the fitted scale-heights to the 
scale-heights derived from observations of edge-on galaxies, as used by DMSvi (EqjJ^. The different lines are for MOND with -7 = 1 and 2 (solid and da.shed black lines respectively). Note all rows u.se the same i-axis range except the galaxy 
UGC 8196. The second column shows the fitted inclination relative to the inclination from the luminous TF relation (VOl). The prior on inclination is the turquoise curve. The third column shows the fitted M/L and the green curve is the confidence 
range of the deduced M/L from DMSvi (using Newtonian gravity with fitted DM halos). 
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Figure 18. As per Fig[T7] 
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Figure 19. As per Fig[T3 
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Figure 20. As per Fig[T7] 
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Figure 21. As per Fig[T3 
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Figure 22. As per Fig[T7] 
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Figure 23. Best fits (black lines) to the vertical velocity dispersions (first and third columns) and rotation curves (second and fourth columns) in MONO with 
7 = 1. The vertical velocity dispersions found using the same best fit M/L and inclination, but the scale-height derived from observations of edge-on galaxies 
(Ea ll6> are given by the red line. The data points are given in black and are derived from the raw measurements of DMSvi, but are corrected for the best fit 
inclination. The error bars on the vertical velocity dispersion data points are the combination of the systematic and random errors and the error bars on the 
rotation curve data points are 10 kms“^ for each point as discussed in E] The dashed turquoise vertical line shows the bulge radius, below which the 
vertical velocity dispersions are unreliable. 
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Figure 24. As per Fig[23] 
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Figure 25. As per Fig[23] 
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Figure 26. As per Fig[23] 




















